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Abstract 

We propose two minimal extensions of Standard Model, both of which can easily 
accommodate the recent puzzling observations about the excess in Rb, the deficit in 
R c and the discrepancy in the low energy and high energy determinations of a s . Each 
model requires three additional heavy vectorial fermions in order to resolve the puzzles. 
The current phenomenological constraints and the new potential phenomena are also 
discussed. 



PACS numbers: 13.38.Dg, 12.15.Ff, 12.90.+b 



Recently it was reported [III by LEP Collaborations that the measured rate of Z — > 66 
is greater than the prediction of Standard Model while that of Z — > cc is smaller. This is 
quite significant given the impressive confirmation of the standard model by other precision 
electroweak tests at the Z° resonance. Given Rb = T(Z — ► bb)/Th a d and i? c = r(Z — > 
cc)/T had , the discrepancies can be summarized as 





Measurement 


SM 


Pull 


Rb 

R c 


0.2219 ±0.0017 
0.1543 ±0.0074 


0.2156 
0.1724 


3.7 
-2.5 



Here SM stands for the standard-model fit with m t = 178 GeV and itlh = 300 GeV, and 
"Pull" is the difference between measurement and fit in units of the measurement error. 

At the same time, the a s problem becomes more acute with improved precision data from 
Z decays. The strong coupling constant a s extracted from high energy measurements at Mz 
seems to be larger than that from low energy measurements, such as deep inelastic scattering 
and lattice calculations |], |3|, ||. The a s (Mz) calculated from the total hadronic width in Z 
decays is 0.125 ±0.0050 f|. On the other hand, low energy measurements all cluster around 
a s (Mz) ~ 0.11. It seems there is a substantial gap between the two. Although more data in 
the future might eliminate these discrepancies, it is possible that this u Rb-R c " plus a s crisis 
are indicating the same new physics beyond the Standard Model. 

Several extensions |5], [j], |8| of the Standard Model have been proposed to address these 
puzzles. In these models, one-loop corrections to the Zbb vertex from the non-standard 
sector will enhance the b quark partial width. With the hadronic total width also enhanced, 
the QCD corrections needed to fit the observed total width is reduced. Thus the observed 
data point to a smaller a s than that in the Standard Model, as favored by low energy 
measurements. However, these attempts all fail to account for the large R c deficit. In 
addition, the first two scenarios might be in potential conflict with top quark decay 0. 

More recently, two papers [0, |ll[] pointed a new direction in the extensions of Standard 
Model which may provide a simple solution to the above discrepancies. Both papers sug- 
gested to resolve the discrepancies by introducing new vectorial fermions that mixes with b 
and/or c quarks. The mixing will reduce or enhance the couplings of the quarks to Z boson 
depending on the gauge quantum numbers of the new fermions. We shall call this class of 
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solutions 'vectorial fermionic solutions' to the puzzles. In Ref.[ll], only a vectorial pair of 
singlet is introduced to reduce the partial width of cc. This could solve the R c puzzle while 
leaving the Rb puzzle only slightly ameliorated. On the other hand, Ref . [|K| , a vectorial pair 
of singlet plus a vectorial pair of triplet are added to resolve both puzzles at the same time 
at tree level. As a price of solving both problems, Ma's model also reduces the prediction for 
the total hadronic width T had and thus renders a surplus in the observed leptonic branching 
ratio Ri = r^/I^, which can not be accommodated by assuming a smaller a s . In this paper 
we propose and analyze two minimal extensions of Standard Model which are nevertheless 
sufficient to resolve the Rb and R c puzzles and simultaneously lower the value of a s ex- 
tracted from Z decay. In the first minimal extension, only a vectorial triplet of fermions are 
needed while; while in the second one, one needs a vectorial singlet plus a vectorial doublet 
of fermions. The first model involves less parameter and therefore has greater predicting 
power. 

We shall start by analyzing the fermion mixing in the general context and then demon- 
strate that our resulting models are indeed the simplest ones of the class. 
In general, the coupling of Z meson with fermions can be written as 

—9—Z»{g f J L ^fL + 9 f R hljR ) (1) 

COS 0\y 

where 

9 f L,R = T^ R -Qf sin 2 6 W . (2) 

The coupling only depends on the weak isospin T 3 and electric charge Q of the fermion. 
Thus mixing with heavy fermions of different weak isospin T 3 could change the coupling of 
quarks with Z meson and the Z decay partial width. Take the partial width into bb as an 
example. Assume that there is a heavy fermion x of charge — | and it mixes with quark b, as 
well as d and s. Denote the mixing matrix among left-handed (right-handed) particles as Ul 
(Ur), which transforms mass eigenstates into gauge eigenstates. We shall specify the weak 
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gauge eigenstates by fields with primes, while those without primes are the mass eigenstates. 



s' 
V 

V x ' ) 



t 





u ds 


u db 


u dx 


u sd 


U ss 


u sb 


U sx 


u bd 


u bs 


u bb 


u bx 


u xd 


u xs 


u xb 


U xx 



\ 



(3) 



V x J 



L^R * ' L^R LfR 

The coupling between mass eigenstate hi and Z° would become 

9l = [Tk\U db \ 2 L + T! L \U sb \ 2 L + T b \\U bb \\ +T* L \U xb \ 2 L -Qsm 2 8 w ] 



(4) 



while g b R equals a similar expression with the subscript L replaced by R. Because the mixing 
matrix U is unitary and quark d, b, s share the same weak isospin T 3 , g b can be written as 

9i, R -T^ R + (T^ R -T^ R )\U xb \l R -Q^n 2 e w . (5) 

The Z partial decay width into bb is proportional to |gfj 2 + \g R \ 2 - 

[T b \ + {T xL — T b 3 L )\U xb \\ -Qsm 2 6 w ] 2 + [T b 3 R + (T 3 R - T bR )\U xb \ 2 R - Qsm 2 6 w } 2 . (6) 

It is different from that in Standard Model. Whether the new fermion will enhance or reduce 
the partial width depends on its weak isospin T 3 . 

Now it is easy to see that we can reduce T c5 by adding a left-handed singlet of T 3 = 
that mixes with Ck|]IU, p"Tf . To increase T bb , a T 3 = —1 left-handed fermion can be introduced 
to enhance \g\\ 2 - A less obvious way is to mix &r, which is of T 3 = 0, with a heavy right 
handed doublet of T 3 = |. 

Next we shall show two minimal extensions of Standard Model in which vectorial fermions 
with the above properties are introduced to resolve both the R b and R c puzzles simultane- 
ously These are the simplest models to accomplish that, with the smallest number of new 
particles, three species of vectorial fermions in both cases. We consider only adding vecto- 
rial fermions since anomalies are cancelled automatically and these fermions could be heavy 
naturally. 

In the first model, only one vectorial triplet is needed. The T 3 = component will reduce 
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T c5 and the T 3 = — 1 component will enhance T b i. The triplet Y can be written as 



Yi 



L,R 



( 5/3 \ 

Vi 

2/3 

y 2 

-1/3 

\y 3 J 



L,R 



with a gauge invariant mass term M Y Y L Y R . The mixing is induced by Yukawa couplings 
between the triplet Yr and left-handed quark doublets. 

UviRt' L <t> + + VuiVl* + + y 3R b' L cf> } , 

UyiRu' L <P + + y2R{u' L <f + d' L <P + )/V2 + y 3R d' L <f>°] . (7) 

In addition we have the ordinary Yukawa couplings in the Standard Model. 

y 3 mixes with the down quarks d, s, b. We'll use the biunitary transformation to diago- 
nalize the 3x3 mass matrix between d, s, b. The mass matrix between D' L = (d' L , s' L , b' L , y' 3L ) 
and D' R = (d' R , s' R , b' R , y' 3R ) then become 

\ 



(d' L ,s' L ,b' L ,y' 3L ) 



( m d 









m s 





£ lV 
£ 2 v 

M Y J 



( 



d' 



R 



\ 



b\ 



R 



D' L M d D' R . 



(8) 



V y 3R J 



M d can be written as 



Mh. 



Md J 



\ 







(9) 



Y J 



with M rf a 3 x 3 matrix, which is diagonal here and J is a 3 x 1 column. It is natural to 
assume that the gauge invariant M Y is much larger than all the other elements of the matrix. 
The diagonalization then takes a simple form. The mixing matrix Ul {Ur) is the matrix 
that diagonalize M d M\ (M ] d M d ). M ] d M d has only one large element at the lower right corner 
with all the other element suppressed by (m/My) 2 . Thus the mixing of y 3R with d R , s R , b R , 
U^i = d,s,b), is also suppressed by {^v/My) 2 and negligible. This is a result of the fact 
that we cannot construct mixing Yukawa couplings between the triplet y 3 L, the singlet q R 
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and the doublet Higgs boson. The mixing between y 3L and b^^s^^d^ is more important. 
Write U L as 

' K R 
S T 

K, R, S are respectively a 3 x 3 matrix, an 1 x 3 column and a 3 x 1 row and T a number. The 



Ur 



(10) 



various elements can be solved in the large My approximation |T2[. In this approximation, 
T is equal to one. K equals the unitary matrix that diagonalizes M^M], which is just unity 
matrix in this case. The column R and row S can also be calculated 



R 



S 



1 



J ] K 



1 



11) 



L Y My My 

The mixing of dn with y 3 L, ie. Ri, Si is approximately ^v/My. As a result, T b i is proportional 
to 

'-^-^ 3 | 2 + ^sin 2 ^) 2 +Qsin 2 ^) 2 , (12) 



with 

To fit the observed Rb, we need 

\s 3 \ 2 



r 



0.0127 ±0.0034 



(13) 



Y , 



(14) 



The charge | quarks will mix with y 2 with the mass matrix 



v My 



(15) 



Note that J is identical to the same column in M^. However M u is no longer diagonal in 
this basis. M u Ml is diagonalized by the KM matrix Vkm- Denote the mixing matrix that 
diagonalizes M u Ml as U' L 

' K' R' 
K S' My 



(16) 



In the heavy My approximation, 



K' ss Vf 



KM 5 



(17) 
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and 



S' = 



My 



•j y w 



KM 



The partial width T c5 is proportional to 



- - -\S: 



I |2 



-sm 2 e w ) + --sin 2 6» 



(18) 



(19) 



,2 2' Zl 3 J V 3 

It is smaller than the corresponding value when there is no mixing. To fit the data, we need 



\S' 2 \ 2 = 0.045 ±0.019 . 



(20) 



Note that S' is related to S through the KM matrix V KM - There is no separate parameter 
for the charge | quarks. The mixing is totally fixed by three parameters £1,2,3 and the heavy 
fermion mass My. 

In the second model, we shall introduce a vectorial doublet and a vectorial singlet. The 
singlet, call it x, has charge | and will reduce R c just like 1/2 in the first model. Choose xr 
to be the only right-handed fermion with the gauge invariant mass term M x xlXr with xl- 
The mass mixing is induced by Yukawa couplings (,iqLi<fixR for % = 1,2, 3. The analysis is 
the same as y<i in the first model. 

The doublet will have the weak hypercharge — f , 



1/3 
-4/3 



L,R 



with a gauge invariant mass term M^ L ^/ R . The Yukawa coupling between \l/ and ordinary 
quarks are 

^0 b R + ^ L 4> s r + £^L0 d R . (21) 

The coefficients £"s need not be the same as the £'s for the singlet x. Therefore more 
parameters are involved in the second model. The mass matrix is 



(d' L ,s' L ,b' L ,*\ 



IL) 



m d 

m s 

m b 

V 6l« &v &v Mm J 



( 



d' 



\ 



R 



D' L M d D' R 



(22) 



Contrary to the previous case, the mixing between ty L with c?l,6l,sl is suppressed by 
(m/M^) 2 . The reason is that we cannot construct mixing Yukawa couplings among the 
three doublets ^[ R , Ql an d Higgs boson. Now the mixing of ^/r with cIr, bR, sr is of the 
order ^v/M^. T b i is proportional to 

2 

(23) 



1 1 

"' "'iw 



+ ^ sin 2 9 W ) + 



3 



To fit the data, we need 



2 



+ - sin &w 



2 V 



0.059 ±0.016. (24) 



Generally, by assuming no non-standard Higgs boson in the theory, if the vectorial 
fermions are triplets or singlets, the effects on gi will dominate because in such case Ur 
is much smaller than Ul. The singlet with the 6-quark charge will only reduce g L which 
is in the wrong direction, while the singlet with the c-quark charge will reduce g c L as data 
require. In the case of a triplet, it can either increase or reduce both g L or g c L depending on 
the hypercharge. 

On the other hand, if the new vectorial fermions are doublet the effects on gR will 
dominate while those on gi are largely unchanged. To increase g b R we need the new vectorial 
"down- type" quark to have T 3 = +|. To reduce g c R we also need the new vectorial "up-type" 
quark to have T 3 = +~. 

From these arguments, it is straightforward to show that the two models we have are 
the ones with a minimal number (three) of new vectorial fermions. If one allows four new 
vectorial fermions, there are also two interesting models that can be considered. One of them 



is adopted by Ma in Ref.|T0[ and the other model uses two vectorial doublets: one doublet 
with Y = -§ to increase g\ and another doublet with Y — | to reduce g c L . We shall not 
discuss these non-minimal models in details. 

In these vectorial fermion models, tree-level flavor-changing-neutral currents (FCNC) 
will in general arise since quarks mix with fermions of different weak isospin. Next we shall 
analyze the FCNC constraints, especially from the kaon decays. 

Because of GIM mechanism, there will be no FCNC if the heavy vectorial fermions have 
the same weak isospin T3 as the quarks they mix with. In the first model, the only component 
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in the neutral current that will generate tree level FCNC in the Kaon decay is — ^'zlI^Uzl- 
It will give rise to a FCNC vertex involving mass eigenstates d and s. 

- 1(^)1^)^^ , (25) 



with 



{Ul)ai — Sx ~ , (U L ) 42 = S 2 ~ ~ S 2 . (26) 



Here ^ 2 f/M is fixed by to be about 0.2 from Eq.(pO|). Thus the coefficient of the FCNC 
vertex 

- -^-^h-fsj. (27) 
2 cos Uw 

is of the order 0.2 x ^v/My- The kaon decay K~l —>■ restricts this coefficient to be 

< 3.1 x 10- 5 [|13[. Take M Y ~ 200GeV as an illustration. The bound for is < 32MeV. 
Given the d quark mass of about 10 MeV, the constraint is still quite natural. If the vectorial 
fermion is heavier, the constrain will be even looser. In the second model, there are more 
parameters involved. The Kaon FCNC constraints will now impose a limit on £i£2- But now 
£2 is no longer fixed by R c fitting, which is related to £ 2 instead. 

In our models, the strong coupling constant extracted from Ri is different from Standard 
Model. While T b i is enhanced and T c5 reduced as the experiment indicates, T uU , T dd - and 
r sg change as well. In the first model, the changes in T uii and T dd are small. Their mixing 
with Y is determined mainly by £1, which is constrained by Kaon FCNC limit. However the 
enhancement in T sS is quite sizable. It is dictated by the necessary change in T c5 to fit R c . 
Overall, without QCD corrections is enhanced. Thus the a s [Mz) extracted from Ri 
using our first model is smaller than using Standard Model since a smaller a s gives smaller 



QCD enhancement corrections^]. Note that Rf, and R c are insensitive to a s . In our model, 
the a s (Mz) calculated from Ri is 0.098 ± 0.007 ± 0.005, with the first error coming from the 
uncertainty in the mixing. It is now consistent with the low energy measurements. In the 
second model, T s is enhanced just like in the first one. Thus the a s extracted from Ri will 
be smaller than 0.125. 

In contrast, Ma's model omits the mixing between the heavy fermion and the s quark. 
Thus T had is reduced since the absolute deviation of R c is larger than that of Rb- Extracted 



from a smaller prediction for Ri, the strong coupling constant becomes 0.18[T3], even higher 



than the original high value of 0.125. This led Ma in his paper [10 to assign the heavy 
fermion a relatively small mass M x < 72GeV so as to open a new channel for the Z boson 
decaying into this heavy fermion. In our models, the a s puzzle is resolved because of the 
enhancement in R s . Experimentally it may be a challenge to measure this R s effect. If this 
could be done, it will be the most direct test of our model. 

The forward and backward asymmetries A b FB and A C FB will also be affected by the mixing. 
The prediction in the first model agrees well with the experimental measurements, as shown 
by Ma[[K]]. In the second model, the asymmetry A b FB is different from the first one. It is 
equal to 0.0980, which agrees with the observed value 0.0997 ± 0.0031. 

One may wonder how these new vectorial fermions can be accommodated in a Grand 
Unified Theory. The vectorial triplet Y can be found is (3, 1, 15) of SU(2) L x SU(2) R x S77(4) 
which in turns can be found in 210 multiplet of 5O(10). The vectorial doublet \I/ is a bit 
harder to accommodate. It can be found in (2, 1, 20) + (2, 1, 20) SU{2) L x SU{2) R x SU(A) 



which in turns can be found in 144 + 144 of 5O(10) or 650 of E e . 
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